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Light-front Hamiltonians for heavy quarks and gluons 

Stanislaw D. Glazek ^ * 

^ Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, Poland. 

A boost-invariant light-front Hamiltonian formulation of canonical quantum chromodynamics provides a heuris- 
tic picture of the binding mechanism for effective heavy quarks and gluons. 
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1. INTRODUCTION 

Light front (LF) Hamiltonian operators for ef- 
fective particles in quantum field theory (QFT) 
provide a new path toward understanding of 
hadronic structure and interactions. The origin of 
new thrust is a far reaching simplification of the 
approach in comparison to the standard form of 
dynamics: boosts are kinematical, instantaneous 
potentials can exist on the front hyperplane with- 
out contradicting the assumption that no physical 
effect can propagate faster than light, and one can 
avoid the problem of the vacuum structure. Ro- 
tational symmetry is dynamical and poses prob- 
lems, but we have a lot of intuition about ro- 
tational symmetry. It seems simpler to obtain 
rotational symmetry in the LF Hamiltonian ap- 
proach than to solve the problem of dynamical 
boost symmetry or the vacuum problem in the 
standard approach. The price one has to pay 
for these gains in the LF Hamiltonian approach 
is a heavy-duty renormalization group procedure 
and a scheme to finesse a leading approximation 
around which one can develop a calculation of cor- 
rections. I describe a heuristic picture of binding 
for heavy quarks and gluons that begins to emerge 
from application of these tools in QCD. I will fo- 
cus on results for gluonium and heavy quarkonia. 

The binding of quarks and gluons occurs typi- 
cally above threshold, which means that the sum 
of masses of the constituents is smaller than the 
mass of the bound state. How can such effect 
occur in a relativistic quantum theory? It does 
not happen in QED. So, how can it happen? A 
boost-invariant light-front Hamiltonian formula- 
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tion of QCD provides a heuristic picture in which 
a binding between heavy effective quarks and glu- 
ons may arise above threshold because the eigen- 
value condition for a renormalized Hamiltonian 
operator may include a large positive contribu- 
tion from mass counterterms that can be compen- 
sated by exchange of gluons only in color-singlet 
states of limited size. One can conceive a sub- 
tle calculation without ever worrying about any 
non-trivial structure in the vacuum. But there 
is a trick needed to finesse effective interactions: 
one introduces an ansatz for a gluon mass gap 
in the Fock states that contain gluons in addi- 
tion to the dominant sectors. Fortunately, the 
effective theory is not very sensitive to the ansatz 
and one can propose a scheme of successive ap- 
proximations that may in principle replace the 
gap ansatz order by order with true interactions, 
"true" meaning implied by the relevant theory. 
It is encouraging that even the first approxima- 
tion exhibits a considerable degree of rotational 
symmetry. The key reference for this talk is 

2. KEY POINTS 

The framework for building a constituent pic- 
ture for hadrons using LF Hamiltonians has sev- 
eral ingredients. The whole procedure begins 
with a canonical Lagrangian of QCD, but the ac- 
tual Hamiltonian operators one works with are 
very different from the bare canonical ones. 

1. In the momentum range where the bind- 
ing mechanism works, the coupling constant a = 
g^/^Air) is comparable to 1. By this I mean 
that the relevant Hamiltonian contains a coupling 
constant g in its vertices, but the constant is a 
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function of a renormalization group parameter A. 
When A ^ oo, the constant g tends to (it is 
never for finite A). In fact, the entire Hamilto- 
nian as an operator is a function of A, 

H = Hx. (1) 

The point is that one needs a method to find 
Hx. In QED, one can try to deduce a Hamil- 
tonian from ^-matrix theory, guessing potentials 
from scattering amphtudes within a perturbative 
expansion for very smaU coupHng constant. In 
QCD, one expects confinement and if this ex- 
pectation is taken seriously into account, includ- 
ing the large value of a in the binding mecha- 
nism, there is no precise link a la QED between 
the S'-matrix for hadrons and quark and gluon 
forces responsible for the binding phenomenon. 
In order to develop LF Hamiltonians for quarks 
and gluons, I will use a renormalization group 
procedure for effective particles (RGPEP). The 
RGPEP provides an expression for Hx order by 
order in an expansion in a small coupling with- 
out reference to the S'-matrix for quarks and glu- 
ons. It also offers a possibility to extrapolate the 
small-coupling results to large values of the cou- 
pling constant because the extrapolation is for 
the Hamiltonian operator, not for the observables 
(see points 5 and 6 below). But the Hamiltonian 
contains finite parts of the ultraviolet (UV) coun- 
terterms that are unknown, and the only possi- 
bility to find those unknown finite parts is to con- 
sider observables for hadrons, including symme- 
tries such as rotational symmetry in decay ampli- 
tudes. 

2. The Hamiltonians Hx must be clearly re- 
lated to the canonical theory. The RGPEP equa- 
tions for evaluating Hx in perturbation theory 
(see e.g. PP) are based on a unitary change of 
basis in the space of operators, 

qx = UxqcanUl, (2) 
Hx{q\) = [Hcan + HcT]reg (qcan) ■ (3) 

This is a similar idea to Gell-Mann's current- 
constituent relationship \'2l^ studied by Melosh 
0], except that now it is built in a dynamical 
scheme of QCD using RGPEP 0. The symbol 
qx stands for creation or annihilation operators 



for effective quarks and gluons, and qcan stands 
for canonical operators. Hcan denotes the canon- 
ical LF QCD Hamiltonian with a regularization. 
An admissible regularization must be imposed on 
the transverse and longitudinal relative momenta 
of interacting particles, those that are invariant 
with respect to the Poincare transformations that 
do not take four-vectors out of the LF hyperplane 
= 0. HcT denotes the required counterterms. 
The RGPEP determines how to find the struc- 
ture of HcT order by order in perturbation the- 
ory. The regularization and renormalization are 
ab initio in the Hamiltonian, defined en block for 
the operators in the Fock space that can be built 
using either Qcan or qx- 

3. In principle, a Fock-space decomposition 
of physical states contains wave functions that 
extend up to the cutoffs introduced by regular- 
ization. Can a Hamiltonian formalism in a Fock 
space with not explicitly covariant regularization 
lead to covariant results? When Dirac introduced 
the front form of dynamics 0, he reduced the 
problem to finding 10 generators that satisfy the 
Poincare algebra: 

[P",?"] = 0, (4) 

[pf'^ArP] = ligf^'pP -g'^Pp-'), (5) 

[M^^M'""] = iigPPM'"' - gP'^MP" 

+ g^PMP" - g'^'M^'P) . (6) 

Can one seek such 10 generators in QFT using 
RGPEP? The seven (one more than in the stan- 
dard approach) kinematical generators: P^, 
M+-L, and J^, do not require regulariza- 
tion and are the same as in a non-interacting the- 
ory. But the three dynamical generators: H = 
P^ and M""*", require regularization and renor- 
malization. So far, it has only been demonstrated 
in a scalar theory up to second order in the cou- 
pling constant [J] that RGPEP produces a re- 
quired solution. In that case one obtains com- 
mutation relations of the form (g is the charge) 

[JlJi] = ie,,kJt + o{g^). (7) 

There is no reason to expect that the method can- 
not work to higher orders. Dirac has observed 
that the challenge of an interacting theory is to 
unfold the commutators that contain products of 
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interactions. Such products appear in second- 
order expressions and a solution exists. Fig. ^ 
illustrates the structure of the first- and second- 
order terms in the creation operators for effec- 
tive particles. The zeroth-order term is a^an itself. 
The Dirac problem in QCD is much more complex 




qcam and (2) an effective Hamiltonian Hx and 
effective operators qx. In each order of pertur- 
bation theory, the result for the S'-matrix is the 
same, provided that the relationship between qcan 
and qx, and between initial H and Hx, is calcu- 
lated up to this order. 

It is important that Hx for effective particles 
contains form factors in interaction terms. Thus, 
it resembles a non-local theory of the type known 
in relativistic nuclear physics, where nucleons in- 
teract through exchange of mesons. In distinc- 
tion from the nuclear models, however, Hx origi- 
nates from QFT and the form factors result from 
a renormalization group procedure. An example 
of old-fashioned diagrams that contribute to an 
amplitude of the type e~^e^ — > hadrons is shown 
in Fig. 13 Note that the integrals over momenta 



Figure I. Structure of the first- and second-order 
terms in in a scalar theory of Ref. 0. The thick 
lines on the left-hand sides of the vertices denote 
canonical creation operators, and on the right-hand 
sides canonical annihilation operators. The thin lines 
indicate how the momentum of an effective particle 
enters the structure. 



than in a scalar model theory. But it cannot be 
addressed in the RGPEP procedure without bet- 
ter understanding of the effective dynamics de- 
scribed by HxQCD because of additional small-x 
singularities in QCD that are not under control 
of the UV renormalization group procedure (see 
below) . 

4. A method for deriving Hamiltonians Hx 
must produce operators that in the case of ex- 
tremely small coupling constant must deliver a 
covariant scattering matrix in perturbation the- 
ory. Wi§ckowski settled the issue in a I- 
loop example of cj)^ theory in 5-1-1 dimensions, 
which is asymptotically free (AF) in perturba- 
tive sense. His more general theorem states: The 
same ^'-matrix for scattering of physical parti- 
cles can be obtained using (1) a bare Hamiltonian 
[Hcan + HcT]reg^ ^^'^ representing the in and out 
particles with creation and annihilation operators 



Figure 2. Counterterms found using RGPEP pro- 
duce a covariant amplitude of the type e+e~ — > 
hadrons, Spa for transition from the incoming state 
a) to the outgoing state |/3), in an old-fashioned one- 
loop calculation with two wave function renormaliza- 
tion factors 1/ in front 0. 



in the loops are limited in a non-covariant way be- 
cause the momenta are already limited once and 
for all in a non-covariant way in the Hamiltonian. 
Nevertheless, the result for the scattering ampli- 
tude is the same as obtained using Feynman dia- 
grams because the counterterms were found from 
RGPEP and all they do is to remove the effects 
of regularization [5] . This is interesting also from 
the point of view of the world-sheet program de- 
veloped by Thorn for planar diagrams |9ll0j . 

5. How could one calculate bound state masses 
using This can be done like in atomic 

physics, where the Coulomb potential is only of 
formal order and still describes a giant variety 
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of bound states. As soon as RGPEP produces 
Hx{gx), one can tackle the eigenvalue problem for 
Hx. A special feature of RGPEP matters here: it 
works in perturbation theory without ever cre- 
ating small energy denominators (this feature is 
built in the design of RGPEP following the prin- 
ciples formulated earlier in |11I12) ). Therefore, 
as long as A is above the size of momenta that 
matter in the binding mechanism, the resulting 
Hamiltonians can provide relatively small matri- 
ces, called windows, whose diagonalization pro- 
duces eigenvalues for the bound states. But the 
closer one approaches the scale of binding, the 
higher order calculation is required and the larger 
is the required coupling constant gx. So, the win- 
dow must be of the right size for the procedure 
to work. This aspect of the RGPEP scheme was 
carefully tested numerically in an exactly soluble 
model with AF A benchmark method of al- 
tered Wegner's equation 14 15 16 produces Fig. 
|31 which shows how successive orders improve the 
accuracy of evaluation of the window in terms of 
its bound-state eigenvalue. Note that the param- 
eter A must be near the scale of binding for a small 
window to work. The model also indicates that a 
low order perturbative window, obtained for arbi- 
trarily small coupling (or arbitrarily small Aqcd 
in the RGPEP scheme) can be extrapolated to 
realistic values of gx and then the window ren- 
ders a good approximation for the bound-state 
eigenvalue of the whole theory. In future, after 
the boost-invariant RGPEP leads to identifica- 
tion of dominant terms via perturbation theory, 
the Wegner equation, or an altered version of it, 
can be applied to calculate window Hamiltonians 
for large coupling constants in specific cases with- 
out using perturbation theory. 

6. In the case of QCD, in the effective-particle 
basis in the Fock space, a gluonium state can be 
written as 

= \gx9\) + \9\9X9\) + (8) 
and a heavy quarkonium as 

= \QxQx) + \QxQx9x) + -- (9) 

But why should such expansion converge at all, 
especially when the coupling constant is compa- 
rable with 1? The reason is that the growth of 



gx when A is lowered using RGPEP can be com- 
pensated by the narrowing of form factors fx in 
the interaction vertices. In other words, the in- 
teractions in a Hamiltonian with a small A die out 
so quickly for large energy changes that they are 
effectively not strong enough to spread probabil- 
ity to sectors with many effective particles, even 
when the coupling constant itself is large |,17j . 
The form factors of RGPEP solve the problem 
that the interactions among effective constituents 
must be strong and at the same time the number 
of the effective constituents must be small (as in- 
dicated by the success of the constituent quark 
model). The situation is akin to nuclear physics 
with practically fixed number of nucleons. 




Figure 3. A bound-state calculation in an asymptot- 
ically free matrix model: a small window Hamiltonian 
matrix is evaluated in six successive orders of pertur- 
bation theory in RGPEP (using altered Wegner's flow 
equation) and then diagonalized non-perturbatively 
using a computer. The exact bound-state eigenvalue 

is -1 [13. 



7. The next key point is that in order to solve 
the eigenvalue problem for Hx in QCD, one has to 
truncate it in the number of effective constituents 
anyway. But in order to do it in a systematically 
improvable way, I assume that there exists a shift 
in the energy of gluons due to non-abelian poten- 
tials that come out from RGPEP in QCD (there 
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are no such potentials coming out in QED, cf. 
[TH|) and that the shift can be approximated in 
the form of a mass gap function for gluons. A 
well-known operation R "19' produces then an ef- 
fective Hamiltonian in the dominant sectors: two 
effective gluons in a gluonium, and a pair of ef- 
fective quarks in a heavy quarkonium. Let me 
follow Maslowski's analysis of gluonium [2111 , as 
an example of the application of RGPEP (the ex- 
ample is analogous to the quarkonium case de- 
scribed in 1 , although gluons demand a more 
advanced analysis because one cannot employ 
a non-relativistic approximation to begin with). 
Matrix elements of the effective Hamiltonian be- 
tween states of two effective gluons \g\g\) labeled 
1 and 2 with relative momentum k, are of the form 



Hkk' = {Ek + EcTk)Skk' 

1 



Wkk' 

1 



Eko — E, 



Ek'o — E, 



(10) 

Yqk', 



where Ek is an eigenvalue of kinetic energy oper- 
ator Ta, equal (fc^^ -t- nil) / {xiX2P'^) , '^^ the 
effective mass of gluons in T\, Eqt is the contri- 
bution of the mass counterterm for gluons, Wkk' 
is the effective interaction generated by RGPEP, 
the sum stands for summation over the effective 
three-gluon basis states that are coupled to the 
two-gluon basis states by the interaction term Y 
in H\, and Eq^ is the LF energy of the three 
effective gluons, each of which has a mass gap 
ansatz in place of m^. The mass gap depends 
on the relative momenta of 3 gluons and can be 
described using a Hamiltonian term 



^.-^E/[123]E^^|123)(123| 



(11) 



where [123] means the LF measure in the integra- 
tion over momenta of the three gluons and the 
sum extends over their polarizations (only two 
transverse ones) . 

8. The trick with the mass gap ansatz is to 
make it correctable order-by-order according to 
the following rule One can write H\ ^ T+V 
and add and subtract the gap term T^, changing 
nothing. But one is free to do it introducing a 
ratio of the coupling constant to its physically 



correct value for a A that one finds suitable to 
work with in this scheme: 

2 



Hx = T + T. 



V - 



(^physical 



T 



.(12) 



For a\ — aphysicai, nothing is changed. But in 
a perturbative expansion for small a\ , the added 
term is large while the subtracted term is neg- 
ligible. This is how one can attempt to incor- 
porate the non-perturbative effects in the three- 
body sector and couplings to sectors with more ef- 
fective gluons in the first approximation. The ef- 
fect of in second-order calculation of Hkk' will 
be replaced by the actual interactions in fourth- 
order RGPEP calculation plus new small ansatz 
correction that will be correctable in higher or- 
der, and so on. The problem one solves this way 
is that unless the sector with 3 gluons is sepa- 
rated by a gap from the sector with two gluons, 
it is not legitimate to use perturbation theory in 
order to account for the coupling between the sec- 
tors. Most probably, in higher order analysis, the 
ansatz for a gap will be pushed away to sectors 
with more effective gluons. Nothing can be said 
for certain yet because the 4th order calculation 
has not been completed. Note, however, that if 
the mass gap ansatz approximates the true QCD 
interactions well, there will be only a little change 
in the leading, approximate picture associated 
with transition from the case with the ansatz term 
Tfj^ and without the interaction V — T^ to the case 
with V and without T^. But in both cases one 
can consider extrapolation of the window Hamil- 
tonians in expansion in powers of the coupling a\ 
from arbitrarily small values to aphysicai that is 
comparable with 1, since the RGPEP form fac- 
tors prevent the interactions from blowing up at 
large momenta for large couplings. 

9. How does the binding mechanism work? 
The mechanism I am talking about can be consid- 
ered an extension of the seminal work of Lepage 
and Brodsky [221 to the domain of small momen- 
tum transfers. They were working on exclusive 
processes with a large momentum transfer us- 
ing diagrammatic rules for scattering amplitudes, 
while the RGPEP allows me to consider what 
happens in the Hamiltonian eigenvalue equation 
(with a gap ansatz) when the momentum trans- 
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Figure 4. The mechanism of binding above thresh- 
old in gluonium. The upper line refers to how the 
mass counterterm for gluons is fixed using a single- 
gluon eigenvalue equation for Hx, and the lower line 
refers to what happens in the eigenvalue equation for 
globally colorless states of two effective gluons, see the 
text. The gluon lines are fat because they represent 
effective gluons. 



fer is in the range of the binding mechanism. 
Of course, Lepage and Brodsky considered usual 
mesons and I discuss here gluonium, but this is 
not essential and I will discuss a case of quarko- 
nium later. Let me explain how the binding above 
threshold emerges using Fig. 01 This is the mech- 
anism that Maslowski used in his calculation [201 ■ 
The main difference in comparison with the ear- 
lier work of Allen and Perry IT , besides the ex- 
plicitly perturbative RGPEP which is different 
from coupling coherence employed by Allen and 
Perry, is that RGPEP provides a boost-invariant 
eigenvalue equation for the mass of a gluonium in 
arbitrary motion, the basis states are built from 
effective gluons instead of the canonical ones, the 
form factors strongly limit changes of invariant 
masses (instead of the changes of that are 
associated in with interaction terms that de- 
pend on spectators), and one has to consider the 
coupling to three-gluon states in order to cancel 
small- a; divergences (3-gluon states could be ne- 
glected entirely in pSjl. 

First one considers an eigenvalue equation for 



a second-order Hx for states with quantum num- 
bers of a single gluon, see the first line in Fig. 
01 One assumes that the mass counterterm is 
such that if this eigenvalue equation were solved 
in perturbation theory to second order in gx (first 
order in ax) the physical gluon mass would come 
out equal 0. This condition produces an effec- 
tive gluon mass which is UV finite but contains a 
positive divergence due to small- a; singularities (a 
blob on the single-gluon line marked with letter g 
in Fig. Q . The divergence is regulated in the ini- 
tial LF QCD Hamiltonian. When one considers 
the same single-gluon eigenvalue problem beyond 
perturbation theory, and inserts a mass-gap in 
the two-gluon sector (now colored, different from 
the colorless states considered in the gluonium 
case below), the negative self- interaction (second 
term in the first line marked with g in Fig. 
is not able to work as strongly as in the case 
of perturbative, massless gluons, and the small- 
X divergence in the effective gluon mass (coming 
from the counterterm) is not compensated: the 
gap blocks the gluons from compensating it be- 
cause they cannot have the required small x easily 
when they have a mass that vanishes too slowly 
for small x. As a result, a single gluon eigen- 
state may have an infinite mass in the limit of 
removing the small-x regularization (Smj^ in Fig. 
E diverges) . 

Next one considers the eigenvalue equation for 
colorless states of two effective gluons and ob- 
serves that the gluon mass gap in a colorless state 
may vanish much faster for small x than in states 
with color. Assuming that it is so, one can ob- 
tain a cancellation of the small-x divergence and 
only a large self-interaction effect is left (finite 
Smj^ in the lower line marked gg in Fig. 0)), sen- 
sitive to the behavior of the mass gap ansatz at 
small X. But there is also an exchange term which 
contributes also a finite but large negative inter- 
action when X of the exchanged gluon is small. 
This exchange term can compensate the large 
self-interaction effect, but not entirely. There 
is a kind of quadratic potential well developing 
around a large minimum whose scale is related to 
the values of A and ax (marked as osc. in Fig. 
21), in addition to a Coulomb force. Of course, 
there are complex spin factors involved and it is 
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a highly non-trivial task to solve the equations 
numerically, see [201 ■ The bottom line is that the 
renormalized self-interactions along the LF in Hx 
build up the gluonium mass high above thresh- 
old of 0. The mechanism in heavy quarkonia is 
analogous but the analysis can be pushed analyt- 
ically much farther because one can exploit the 
non-relativistic limit when the quark masses are 
much larger than A. This mechanism can be con- 
sidered a generalization of the 1+1 dimensional 
model of 't Hooft _24j to 3-1-1 dimensional QCD. 

3. GLUONIUM 

In order to illustrate what comes out from LF 
Hamiltonians for effective gluons in pure gluody- 
namics, let me draw on Ref. Maslowski [201 
considered a class of the ansatz masses for every 
gluon in the three-gluon sector, i = 1,2,3, given 
by the same formula 

^2 ^ Xi{xiX2X-iY , (13) 

where b ^ \ ^ 2 GeV, and 6 ~ 0.2 are constants. 
A typical dependence of the smallest gluonium 
mass on a is shown in Fig. [S| When the cou- 
pling constant increases, the state collapses. But 
it is known from the beginning that the eigenvalue 
equation with only two or three effective gluons 
cannot be valid for states with masses much larger 
than two times an effective gluon mass. The aver- 
age value of the gap ansatz turns out to be about 
1.5 GeV in all interesting cases. Therefore, it is 
encouraging that there exists a range of couplings 
for which the gluonium mass is large and increases 
with the coupling. 

An example of a spectrum of smallest masses 
from [201 is shown in Fig. [01 and in a table: 



M(j, = 0) 


M{\]A = 1) 


M(|j.|-2) 


1.73 


2.25 


1.97 


2.01 


2.63 


2.54 


2.06 


2.64 


2.6 


2.41 


2.66 


2.67 


2.41 


2.85 


2.83 



These results represent a typical output from 
Hx- The example was obtained using 9 radial 
wave functions and 15 spherical harmonics. The 
results were stable with respect to changes of 



these two numbers. The resulting spectrum is 




Figure 5. Typical dependence of the smallest 
gluonium mass (in GeV) on the coupling constant 
a between 0.15 and 0.5. 



slightly denser than on the lattice but otherwise 
it appears to be qualitatively similar. However, 
there is so far no clear way for associating the 
eigenstates with spin quantum numbers. There 
is no clear degeneracy into multiplets correspond- 
ing to rotational symmetry. Although the low- 
est masses are not widely different, and one may 
think that the violation of the multiplet structure 
is not significant given the crude nature of the cal- 
culation (cf. 23 ) , the question of how to improve 
rotational symmetry is not answered. But it is 
certain that a simple mass gap ansatz of Eq. H13|l 
has no a priori reason to produce true degeneracy 
of the masses. 

Regarding the stability of the results versus 
changes of A one can say that when ax is forced 
to vary with A according to the perturbative for- 
mula with no quarks, the results are not varying 
significantly with A. But one can change the rel- 
ative order of masses of \jz\ — 0, 1,2 eigenstates 
by making about 10% changes of the values of A, 
a, h and 5. The overall conclusion is that the glu- 
onium picture in the boost-invariant LF Hamil- 
tonian approach appears surprisingly reasonable 
even in the crude first approximation. Readers in- 
terested in more details should consult Ref. |20| . 
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Figure 6. Gluonium masses: dashed lines are taken 
from lattice calculations (55]. Bold lines are obtained 
from H\. The mass of the state labeled 0^"'' was fitted 
to the lattice result of 1.73 GeV. The parameters in 
Hx are A = 1.92 GeV, b = 2X, S = 0.2, q, = 0.44. 



4. QUARKONIUM 

Theoretical aspects of the heavy quarkonium 
picture emerging from LF Hamiltonians for QCD 
with one heavy flavor and gluon mass gap ansatz 
are described in Ref. PP| and do not need to be 
repeated here. The main feature is that all de- 
tails of the mass ansatz for effective gluons in 
the sector \Q\Q\gx) disappear from the effective 
Hamiltonian in the \QxQx) sector and the result- 
ing second-order eigenvalue problem turns out to 
be exactly rotationally invariant. The first cor- 
rection to the Coulomb potential is a force re- 
sembling a harmonic oscillator. 

The fact that a mass gap ansatz leads to an 
oscillator-like interaction term which respects ro- 
tational symmetry already in the 2nd order anal- 
ysis using ~ 1 does not seem accidental. This 
result appears almost independently of all de- 
tails of the ansatz because momentum transfer 
squared, q^, in the terms that are sensitive to 
small- a: of the gluons, is limited by the RGPEP 
form factors f\ to so small values that the ratio 
/i^/(<Z^ + M^) is practically 1 for any reasonable 



ansatz. In addition, it seems likely that the same 
result comes out also as a part of the genuine 4th 
order calculation (not completed yet). In the 4th 
order, the ansatz term ~ 1 cancels out for 
large physical a The next term comes from 
QCD interactions order a in the 3-body sector 
and this is how the actual gap may show up. But 

in a system dominated by the Coulomb force 
(the case of quark masses very much larger than 
Aqcd) continues to be formally on the order of 
the strong Bohr momentum squared, or order 
which is much smaller than /i^ when fi^ ^ a. 
If this observation is confirmed in 4th order cal- 
culations, the spherical symmetry of the leading 
oscillator term identified already in the 2nd or- 
der using the ansatz for /i^, may be a necessary 
consequence of the rotational symmetry of QCD. 

One should observe two general arguments that 
support the harmonic result in the first approxi- 
mation. One argument is that the combined effect 
of the self-interactions and the exchange leads to 
binding above threshold which emerges around a 
minimum of the mass squared operator derived 
from the LF Hamiltonian and every func- 
tion around a minimum is in the first approx- 
imation a quadratic one. The question is not 
so much why the potential comes out quadratic 
but rather what the spring (not string) constant 
is. Quite general Coulomb scaling argument im- 
plies that the spring constant is on the order of 
k ~ aX^ /m^. But if one assumes that A ~ sfmri^ 
the resulting harmonic potential scales with a ex- 
actly as the Coulomb eigenvalue problem does. 
This means that such harmonic force may always 
be of a fixed relative magnitude to the Coulomb 
part of the interaction if QCD develops a mass 
gap for the effective gluons. The other argu- 
ment is that when the distance between the heavy 
quarks increases, a linear potential may develop 
between the quarks at the expense of creating 
additional gluons with a mass gap if the poten- 
tial grows faster than linearly |26I27| - it becomes 
energetically more favorable to create a definite 
number of gluons per unit of length than to al- 
low the potential energy to grow without creating 
additional gluons. And the first possible integer 
power of the distance is 2. Thus, the harmonic 
force is probably the simplest one that can lead 
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to a string picture in LF Hamiltonians for QCD. 

The simplicity of the LF Hamiltonian approach 
can now be exhibited using a calculation per- 
formed by a freshman on a personal computer 
|28|. When one neglects all spin effects and then 
disregards parameters that play no role in the 
spinless case, the eigenvalue equation for heavy 
quarkonia that one derives from the LF Hamilto- 
nian reads 



'2m 



4a 
37 



r 



^{r) = E^{r), 



(14) 



where m is the reduced quark mass, equal half of 
the quark mass mq. The eigenvalue E gives the 
quarkonium mass through the formula 



M = 2m, 



,\Jl + E/niq , 



(15) 



which is relativistic. 

Stawikowski found that for k = 0.0431 GeV^, 
2m = mb = 5.0006 GeV, a = 0.6029, aU compati- 
ble with the range of parameters a priori possible 
in LF QCD, the masses of T states IS', 2^, and 
35", could be reproduced with precision of 10~^, 
or ±0.05 MeV. Such accuracy is understandable 
because one has three parameters to fit three 
masses (Stawikowski considered a Coulomb po- 
tential multiplied by the factor e^^'"(l — e'^'') but 
for the values ^ = 10"^ GeV and A = 4 GeV that 
he used, these additional parameters were not 
important). The same parameters give Mip — 
9.9144 GeV (the experimental average of masses 
of IP bottomonium states is 9.8884 GeV), and 
M2P = 10.2507 GeV (the experimental average 
of masses of 2P bottomonium states is 10.2519 
GeV). 

Keeping the same parameters k and a and 
changing the quark mass to nic = 1.6068 GeV, 
one obtains the experimental value of the mass 
of J/^ state. Mis = 3.0969 GeV. The 25" state 
obtains the mass K'hs — 3.6628 GeV (the exper- 
imental value is 3.6861 GeV), and the state IP 
obtains Mip — 3.4991 GeV (the experimental av- 
erage of masses of IP charmonium states is 3.4940 
GeV). It is clear that the harmonic oscillator po- 
tential as an admixture to the Coulomb potential 
is a qualitatively acceptable one from the point 
of view of such fit. An extensive numerical study 
including spin effects is under way 29 . 



5. CONCLUSION 

The second-order RGPEP procedure for deriv- 
ing LF Hamiltonians for heavy quarks and gluons 
produces operators that seem to have a chance 
of providing a qualitatively reasonable first ap- 
proximation of the dynamical picture for masses 
of gluonium and quarkonia. Rotational symme- 
try requires intensive studies, especially in the 
case of gluonium (important gluon degrees of free- 
dom appear also in hybrids and there one faces 
considerable constraints due to rotational sym- 
metry [siQj). The picture of binding mechanism 
that one can hold on to in order to build further 
intuitions is based on the possibility to control 
behavior of self-interactions of quarks and glu- 
ons in the effective Hamiltonians with small RG- 
PEP parameter A. The self-interactions largely 
reduce positive contributions from counterterms, 
but not entirely, and the exchange of gluons can 
lead to bound-state masses above threshold. The 
RGPEP procedure can be applied to many more 
cases than only those discussed here and one can 
refine the picture using higher order perturbation 
theory and a whole arsenal of methods for solv- 
ing quantum Hamiltonian problems. An exten- 
sion of the approach to include light quarks will 
be more complicated than the heavy quarkonium 
case. But the example of gluonium suggests that 
the required fourth order RGPEP studies may 
shed new light on the issue of how to approach 
eigenvalue problems including light quarks. So 
far, no vacuum structure was involved in the en- 
tire calculation. 

The rotationally symmetric effective Hamilto- 
nians described here provide an interesting addi- 
tion to earlier studies based on similarity renor- 
malization group and coupling coherence |31I32[ 
I33I34I35| . But it is not clear yet if the RGPEP 



procedure for LF Hamiltonians will ever match 
the accuracy of methods applicable in QED |36l 
I57| or in lattice-based calculations ■38"39". One 
has to be aware that a constituent picture can 
be fitted to data using a large range of potentials 
|4()l41j . On the other hand, if the hypothesis that 
QCD has an infrared limit cycle 42,4.3^ is correct, 
it may turn out that new Hamiltonian methods 
provide a path to understanding of the limit cycle 
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universality required for systematic solution 
of the theory. 
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